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Abstract 



Due to a rapid change of the entropy density s(T) across the critical 
temperature T c of the QCD phase transition, the pressure P{T) and the energy 
density e(T) above T c generally deviate from their Stefan-Boltzmann values. 
We shall demonstrate this both analytically and numerically for a general class 
of s{T) consistent with thermodynamical constraints and make a qualitative 
comparison of the result with the lattice QCD data. Quantities related to 
ds(T)/dT such as the specific heat and sound velocity are also discussed. 
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Revealing the precise nature of the quark-hadron phase transition is one of the central 
issues in recent lattice and analytical studies of hot QCD (see the reviews Q.) Furthermore, 
the behavior of the bulk quantities such as the entropy density s(T), energy density e(T) 
and pressure P(T) as a function of temperature T is relevant to the formation and evolution 
of the quark-gluon plasma in the ultrarelativistic heavy ion collisions planned at BNL and 
CERN. 

Lattice QCD calculations have been providing interesting data on the plasma properties 
above the critical temperature T c []|. Among them are, however, some unexpected features: 
(a) P/T 4 approaches the Stefan-Boltzmann limit very slowly as T increases, (b) e/T 4 has 
a peak just above T c and approaches an asymptotic value from above as T increase, and 
(c) e — 3P 7^ above T c i.e. a large deviation from the ideal gas behavior is seen. They 
are sometimes identified as indications of the non-perturbative nature of the quark-gluon 
plasma. 

In this letter, under generic assumptions on s(T) consistent with the thermodynamic 
inequality, we examine what is the expected behavior of e(T), P(T) and other quantities 
such as the specific heat Cy(T) and the sound velocity c s (T). We show that most of the 
"unexpected" behaviors measured in the lattice QCD simulations |]-[J can be explained at 
least qualitatively as a simple consequence of the rapid increase of s(T) around T c . 

Our basic observation is that all the above quantities are simply parametrized by the 
entropy density s(T) which depends only on the active degrees of freedom and is free from the 
complexities of the vacuum structure of QCD. This is seen by the thermodynamic relations 
at zero chemical potential such as 

P(T) = j\(t)dt, (1) 

e(T) = Ts(T)-P(T), (2) 

and 

Cv = TW (3) 
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(4) 



In eq.(f]), P(T) is normalized by P(T = 0) = 0, which also leads to e(T = 0) = from 
eq.(|2|). The same normalization is adopted also in lattice calculations. 

Since we are interested only in the gross behavior of the thermodynamic quantities and 
not in the precise order of the phase transition |4[], we make a smooth interpolation of 
s(T) between the the hadronic gas at low T and the quark-gluon plasma at high T. Note, 
however, that one cannot make arbitrary parametrization of s(T) since it is constrained by 
the thermodynamic inequality and the Nernst's theorem |5| 

^>0, *(0) = 0. (5) 

A simplest possible parametrization satisfying eq.(E]) reads 
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s(T) = s h (T)w h (T) + s q (T)w q (T), 



(6) 



where Wh{T) = 1 — w q (T) and w q (T) is an increasing function. In this letter, we take the 
following form for w q (T), but our conclusion is not limited to this specific choice of Wh, q {T): 



n (l + tanh (J— / 
Wq ^ T) = m (l - tanh (^)) + n (l + tanh V ' 



Here Sh(T) = 12(7r 2 /90)T 3 and s 9 (T) = 148(7r 2 /90)T 3 are the entropy densities of massless 
free gas with two flavors in the hadronic phase (pion gas) and the quark- gluon phase, re- 
spectively. The interaction between particles and the quark masses are neglected just for 
simplicity. 2T sets the width of the transition region. For T satisfying \T — T c \ ^> T, s(T) 
approaches asymptotically to Sh(T) below T c and to s q (T) above T c . m and n are introduced 
to consider the asymmetric superposition of the two phases around T c , but we shall take 
m = n in the following. 

The equation of state obtained from eq. (|6|) is more general than that of the bag model 
in the sense that (a) we need not refer to phenomenological parameters such as the bag 
constant, and (b) we can treat not only the strong first order phase transition but also 
the crossover by changing T in a thermodynamically consistent way. When T = 0, our 
equation of state is equivalent to the bag model one with a bag constant 4B = (s 9 (T)/T 3 — 
s,(T)/T 3 )T 4 . 

In Fig.l, s(T) with m = n = 1 is shown for T/T c = 0, 0.05 and 0.25. Lattice measure- 
ments show a rapid variation of s(T) in a narrow region of T (with a width ~ 10 MeV) 
0, which is similar to Fig.l with T/T c = 0.05. Inclusion of the perturbative a s corrections 
above T c and the interactions of hadrons below T c modifies the absolute value of s(T) as 
well as its T dependence, but it does not change our conclusions qualitatively. 

Once s(T) is given, it is straightforward to calculate other quantities from eqs. (|l]J^J^,|) . 
In Fig. 2, e/T 4 , P/T 4 and (e— 3P)/T 4 are shown as a function of temperature for r/T c = 0.05. 
They can be directly compared with the lattice results [1|-|3| at least qualitatively. In fact, 
they look quite similar: 

(i) P/T 4 increases rather slowly above T c both on the lattice and in Fig.2. Since 
P(T) is given as an integral of s(T), it is quite natural to expect such a continuous and slow 
rise. When T = 0, P(T > T c ) has an analytic form 

P(T)/P SB (T) ~ 1 - (T c /T) 4 , (8) 

where Psb(T) denotes the Stefan-Boltzmann value. Eq.(||) tells us that P(T)/Psb(T) = 
50% (90%) for T/T c = 1.2 (1.8) independent of the details of the dynamics. If T is increased, 
P/T 4 has even weaker T dependence. Thus the major part of the deviation of P(T) from 
Psb{T) can be accounted for without introducing non-perturbative interactions of quarks 
and gluons above T c . 

(ii) e/T 4 has a peak just above T c in Fig. 2 (note that e(T) itself is a monotonically 
increasing function of T), which is also seen in lattice simulations Qj. In our case, this peak 
is a simple consequence of the rapid increase of s(T) and the slow rise of P(T) above T c 
(see the definition eq.(§)). The width of the peak is correlated with the slow rise of P/T 4 . 
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One can even prove analytically that there must exist a peak around T c from the following 
relation satisfied by arbitrary weight functions Wh, q {T), 

T 5 = T 2 (s h (T)w' h (T) + s q (T)w' q (T)) 

- f T t (s h (t)w' h (t) + s q (t)w' q (t))dt, (9) 
Jo 

where the prime denotes a derivative with respect to T and we have used a fact 
Sh, q (T) = (constant) xT 3 . Under the conditions that Whif) + w q {t) = 1, < w q (0) < 1, 
w'{t) > and w'Jt — > oo) = o(l/t 4 ), (e/T 4 )' is positive (negative) for low (high) T and 
has a zero near T c . These conditions are all satisfied in our case, and as a result e/T 4 
has a peak. One should also note that, for high enough T, (e/T 4 )' ~ — A(T)/T < with 
A = (e - 3P)/T 4 . 

(iii) A(T) = (e — 3P)/T is so called the "interaction measure" and has a peak near T c 
both on the lattice p]|§ and in Fig.2. Again, the rapid increase of s(T) (i.e. the liberation 
of quarks and gluons) is the reason for this peak, which is also seen from the general formula 

A(T) = jT t (s h (t)w' h (t) + s q (t)w' q (t))dt. (10) 

A vanishes for T = and oo and has a peak around T c independent of the details of the 
dynamics. Also the peak becomes sharper and its height becomes higher as T decreases. 
When T = 0, one can rewrite A as A(T > T c ) = (s q /T 3 — Sft/T 3 )(T C /T) 4 which is equivalent 
to the bag model formula without a s corrections A(T > T c ) = 4P/T 4 ||. Note here that, 
if we normalize the pressure by P(oo) = P SB instead of P(0) = 0, the peak of A(T) does 
not arise and the naive Stefan-Boltzmann law is realized above T c . In this sense, the 1/T 4 
tail of A at high T in Fig.2 can be interpreted as an artifact of the normalization and has 
nothing to do with non-perturbative interactions of quarks and gluons in the quark-gluon 
plasma. 

In Fig. 3, c 2 is shown for r/T c = 0.05. Since the equation of state becomes soft near the 
critical region, the sound velocity slows down. This is the reason why c 2 has a sudden drop 
in the narrow region T c — r<T<T c + r. This sudden change is in contrast to the broad 
peak or slow rise of the quantities in Fig.2. The difference comes from the fact that c 2 is 
related to the derivative of s(T) while e(T) and P(T) are related to the integral of s(T). To 
see the effect of the finite pion mass on c 2 below T c , c 2 using Sh(T) with = 140MeV (and 
T c = 180MeV) is shown in Fig. 3 by the dashed line. The effect of finite on the other 
quantities is small since they are small in any way at low T. The heat capacity Cy has also 
a sharp peak at T c , since one needs to supply a large amount of heat to increase T across 
T c to liberate quark-gluon degrees of freedom. 

Although we need not refer to the vacuum parameters such as the bag constant in our 
approach, it might be instructive to introduce an effective bag constant above T c defined as 
a deviation of the pressure from its Stefan-Boltzmann value, 

B(T > T c ) = -(P(T) - P SB (T))/T\ (11) 

where we have neglected all the a s corrections. B(T)/T^ is shown in Fig. 4. The asymptotic 
value B(oo) is the one usually used in the bag model. In our case, B(oo) depends on how 
one parametrizes s(T): as T is increased, B(oo) also increases. 
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The main conclusion of this letter is that, whenever the entropy density has a rapid 
change near T c , e/T 4 , P/T 4 and A behave as we know from the lattice simulations. Fur- 
thermore, we need not refer to the vacuum condensate or the bag constant to see this fact. 
We believe that our approach based on the parametrization of s(T) provides a transparent 
and thermodynamically consistent way to study the qualitative feature of the bulk plasma 
quantities. Although we have taken a simplest possible form of s(T), inclusion of the quark 
masses, hadronic interactions, a s corrections and chemical potentials is straightforward. 
Possible non-perturbative effects at T c < T < 3T C [||,|| appear as a deviation from the basic 
curves in Fig. 2 and could be included by adjusting the functional form of s(T). However, 
it is not an easy task to identify true non-perturbative effects from the lattice data, since 
Fig. 2 has already similar behavior with the lattice data, and furthermore the finite volume 
effect and perturbative corrections are still large in the current lattice simulations M. 



ACKNOWLEDGMENTS 



We thank Teiji Kunihiro for useful comments and discussions, and Rob Pisarski for warm 
hospitality at BNL where most part of this work was completed. We also thank O. Miyamura 
and M. Gyulassy for private communications who have already realized a part of the present 
results in the context of the bag model. M. A. was supported by the Director, Office of 
Energy Research, Division of Nuclear Physics of the Office of High Energy and Nuclear 
Physics of the U.S. Department of Energy under Contract No.DE-FG02-93ER40764. T. H. 
was supported in part by the Grants-in-Aid of the Japanese Ministry of Education (No. 
06102004). 



4 



REFERENCES 



[1] See e.g. C. DeTar, |hep-ph/9504325| (1995), and K. Rajagopal, |hep-ph/95043T0| (1995), to 

appear in Quark Gluon Plasma 2, edited by R. Hwa (World Scientific, Singapore, 1995) 

and references therein. 
[2] T. Blum, L. Karkkainen, D. Toussaint and S. Gottlieb, Phys. Rev. D51 (1995) 5153. 
[3] E. Laermann, Nucl. Phys. B (Proc. Suppl.) 42 (1995) 120, and F. Karsch, |hep-lat/95030TO| 

(1995) and references therein. 
[4] For massless two flavors, the chiral transition is likely to be weak first order (F. Brown 

et al, Phys. Rev. Lett. 65 (1990) 2491; Y. Iwasaki et al., ^ep-lat/9505Uil (1995)). One 

can mimic it by choosing a more sophisticated form of Wh, q {T), which, however, will not 

change our conclusions qualitatively. 
[5] L. Landau and E. Lifshitz, Statistical Physics, (Pergamon, New York, 1980). 
[6] H. Satz, Phys. Lett. B113 (1982) 245. 
[7] N. H. Christ, Nucl. Phys. A527 (1991) 81c. 

S. Chandrasekharan, Nucl. Phys. B (Proc. Suppl.) 42 (1995) 475. 
[8] V. Koch and G. E. Brown, Nucl. Phys. A560 (1993) 345. 

D. H. Rischke and W. Greiner, Phys. Rep. (1995) in press, and references therein. 



Figure Captions 

Fig.l: s/T 4 as a function of T for three different values of the width parameter T. 

Fig. 2: The solid line, the dashed line and the dash-dotted line correspond to e/T 4 , P/T 4 

and A, respectively. T/T c = 0.05 is chosen. 

Fig. 3: Squared sound velocity c 2 s as a function of T. Solid (dashed) line corresponds to 
= (140) MeV. 

Fig. 4: "Effective" bag constant B(T)/T 4 for T > T c . Solid, dashed and dash-dotted lines 
correspond to T/T c = 0, 0.05 and 0.25, respectively. 
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